Coherent field emission image of graphene predicted with a microscopic theory 
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Electrons in the mono-layer atomic sheet of graphene have a long coherence length of the order 
of micrometers. We will show that this coherence is transmitted into the vacuum via electric 
field assisted electron emission from the graphene edge. The emission current density is given 
analytically. The parity of the carbon pi-electrons leads to an image whose center is dark as a result 
of interference. A dragonfly pattern with a dark body perpendicular to the edge is predicted for the 
armchair edge whose emission current density is vanishing with the mixing angle of the pseudo-spin. 
The interference pattern may be observed up to temperatures of thousand Kelvin as evidence of 
coherent field emission. Moreover, this phenomenon leads to a novel coherent electron line source 
that can produce interference patterns of extended objects with linear sizes comparable to the length 
of the graphene edge. 

PACS numbers: 79.70.+q, 72.80.Vp, 61.05.jp 



Nano-emitters of various shapes and materials have 
been used for cold field electron emission (CFE) with 
the aim of constructing a practical microelectronic vac- 
uum electron source that may be used in flat-panel dis- 
plays, electron microscopes and parallel e-beam lithogra- 
phy systems. The feature most important for such ap- 
plications is the high aspect ratio of nano-tips enhancing 
the apex field to the point where field electron emission 
can be driven by macroscopic electric fields of about ten 
volts per micrometer or even less. Two successful ex- 
amples of this approach are Spindt-type cathodes |l| and 
single atom field emission tips.@, y| In the context of 
the present paper we should note that in the latter the 
emission volume is the terminating atom so that emitted 
electrons are spatially coherent as proven succinctly in 
in-line electron holography Q The coherent field emission 
of carbon-based nano-emitters has been confirmed with 
the perfect interference fringes produced by the electron 
beam from carbon nanotube CFE. [5|, |(| 

Several groups have demonstrated that graphene does 
show promising CFE properties such as a low emission 
threshold field and large emission current density. 0413 
What we will show in this paper is that the unique elec- 
tronic properties of graphene and its monolayer atomic 
structure result in a novel electron line source with co- 
herent field emission that can produce interference pat- 
terns of extended objects with linear sizes comparable to 
the length of the graphene edge. What makes graphene 
unique is that its two-dimensional electron system is 
well described by a superposition of 2p carbon orbitals 
and that its low energy excitations behave like Dirac 
particles. HHlil Due to the perfect atomic structure, 
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electrons in graphene have long coherent lengths (a few 
micrometers at room temperatures [l7j). We will show 
that the emission electrons will carry this coherence of 
the quantum states of graphene into the vacuum. 

The challenge is to explain/predict the coherent as- 
pect of the CFE image. Because of its wave nature 
coherent emission electrons must have been in a collec- 
tive quantum state in the solid and then be transmit- 
ted through the tunneling barrier into a single particle 
state in the vacuum. To address the coherence aspect 
of CFE, one must incorporate the electronic structure 
of a mesoscopic/macroscopic nano-emitter and deal with 
the screening of the macroscopic applied field and of the 
atomic orbitals along the edges so that one is in a posi- 
tion to calculate the tunneling wave function. Needless 
to say this goes beyond the framework of conventional 
CFE theory as the phase of the quantum states is cru- 
cial here. But it can be done for graphene as shown in 
this paper. We will construct the wave function of the 
emitted electron analytically thus extracting important 
information about the collective quantum state. 

A piece of graphene with the armchair edge (AE) along 
the y-axis is shown in Fig (Ha) . If it rotates by an angle 
of 90 degrees on the plane, the emission edge will become 
the zigzag edge (ZE). We adopt a rectangular plaquette 
of the red lattice as the unit cell (see Figflja)). There 
are four carbons in a unit cell, each is labeled by a — A± 
or B± and referred to as the a-carbons. The energy 
bands are plotted in FigQJb). The CFE experiment is 
schematically shown in FigQ](c-f), with (c) the graphene, 
(d) the anode plane, (e) the electric objective lens and 
(f) the back focal plane(BFP), respectively. A uniform 
macroscopic field Fq is applied in the negative direction 
of the a;-axis. 

The basic microscopic process of the CFE may be un- 
derstood as an electron in the graphene being annihilated 
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and then created in the out-going state in the vacuum. 
This is the transition of a two-dimensional massless ex- 
citon to a three-dimensional elementary electron. The 
low energy excitations obey Dirac dispersion relations 
E(k) = ±Suf|k| , with ± for the positive /negative en- 
ergy bands, vf ~ c/330 the Fermi velocity. 1 1 51] and k 
a wave vector relative to the K point (FigOJb)) on the 
x-y plane. We will focus on the positive energy bands 
which has dominate contribution to the CFE, although 
our analysis can be easily extended to including the neg- 
ative energy bands. The eigenstates are Bloch waves 
along the y-direction but not in the ^-direction due to 
the edge. We have adopted the time-Y parity representa- 
tion such that the eigenstates of the half-infinite graphene 
have either even (s = +) or odd (s — — ) parity in the 
combination of the time-reversal and the y-directional re- 
flection transformations. The eigenfunctions can be ex- 
panded in the Wannier-Bloch hybrid orbitals that in the 
x-direction are localized at each carbon and extended in 
the y-direction, 

x , e ik y y 

where &2 = §f is the reciprocal lattice spacing in 
the y-direction, k l y = k y — tt>2 with I an integer, and 

t Pm,ki( x ' z ) i s defined as the /c^-component of the y- 
directional Fourier transformation of the Wannier orbital. 
The Bloch- Wannier orbital with k l y has y-directional ki- 

(hk l ) 2 

netic energy E y = 2 ^ and the energy for the mo- 
tion on the x-z plane is E p = E(k) — E y . This leads 
to an effective vacuum barrier height W p = W — Eg, 
with W ~ 4.32 eV the work function of graphene, [18| 
for tunneling normal to the edge. [H, H(| For hydrogen 
saturating edges, E y of the AE vanishes at the K point 
if I = while it is 2.8eV for the ZE. Therefore emission 
from AE with k y = k y is dominant. However we should 
note that the transverse wave numbers of edge states of 
the ZE are shifted to the vicinity of the T point if the 
edge is saturated via atoms with p-orbitals (oxygen for 
instance) with local energy similar to the 2p orbital of 
carbon. [21| In such a case CFE from edge states of ZE is 
also possible. In the present paper we will focus on the 
band state emission and only consider the I = mode of 
the AE. The subscript I will be dropped latter on. 

The amplitudes C" m (k) are determined by the bound- 
ary condition and the time-Y parity. Notably, the AE 
will mix states with the same bulk energy and a given 
k y . These degenerate states include excitations from 
both Dirac cones {K and K' in FigOJb)) and states 
with k x and — k x . In the final expression for the emis- 
sion current density, we only need the combination Ck = 
|C#(k) + C*f (k)| = sin |jig 71 with 7 the pseudo-spin 
mixing angle (the phase difference between the two sub- 
lattices) given by tan 7 = fc^/fc^.flSjl and M (N) the total 
number of rows (columns). 

Because the relevant electronic properties of graphene 



can be attributed to the 2p orbital of carbon, the 
edge Wannier orbital should have the form <fk(p) — 
N]g_Zy/ge~ TpP with p = [x, z) and g — \J x 2 + z 2 inside 
the region where the atomic potential is dominant (APD 
region). Here the factor z comes from the 2p orbitals, 
TVk is a normalization factor and r p = i yj2mW p . 

The graphene eigenstates are not exactly stable in the 
applied fields. The emission wave should be found by 
solving the Schrodinger equation. The probability am- 
plitude emitted from -0ks for large times is strictly given 
by the path-decomposition formula. [22T| 

^»(r) = jf GT(ry)^£ ^G(r',r")^ s (r") (2) 

where J„ is the surface integral over an arbitrary surface 

E that separates the graphene and the vacuum, and d n is 
the right normal derivative minus the left one on E. The 
G is the retarded fixed energy Green's function and the 
restricted Green's function G r is defined by a path inte- 
gral over energy-fixed paths that stay strictly outside S 
and do not touch it except at the initial point. The physi- 
cal meaning of © is clear. The emission electron must go 
through the separation surface E on which the amplitude 
Gipks at each point acts as a point source. The operator 
~2m^" represents the mean imaginary flux through the 
surface. The Green's function G r propagates the wave 
from E to the screen. Thus the surface integral gives 
the emission wave as the coherent superposition of waves 
from all points on E. We will choose E as a surface of 
constant potential in the vacuum barrier region such that 
only the least action path has significant contribution to 
G r and the semi-classical approximation (JWKB) is ap- 
plicable. We also require that E locates on the boundary 
of the ADP region so E is invariant in the y-direction yet 
the orbital wave function is still feasible on it. Because 
G(r',r") decreases with |r' — r"| exponentially, only the 
edge carbon orbitals with m = 1 in the expansions (TTJ) 
have significant contribution to ([2]). 

To go further, the vacuum barrier potential in the 
vicinity of the edge is needed. It may be written as 
V{p) — W + V a t{p) + Ve{p) with V a t the short-ranged 
atomic potential and V e (p) the electric potential energy. 
First note that the y-directional momentum hk y is con- 
served due to the symmetry. More complicate interac- 
tions would be involved when the electron is very close 
to the atom core. But the exact form of V a t close to 
the atomic core is not crucial if the wave function is 
known in the APD region as we have assumed. A feasible 
model is V a t = — with P the dipole line density along 
the edge.[IH It is easy to obtain V e outside the APD 
region (referred to as the off- APD region) via the con- 
formal transformation, z + ix — y[z + i(x + h)) 2 + h 2 
with h the height of the graphene. On the p — (x,z) 
plane, V e = —eFox. With the screen far away from the 
graphene, the observed image is produced by the near 
axis (x or 5-axis) electrons in the off-APD region. In 
the near axis approximation (NAA), neglecting terms of 
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order |y , one gets V at ~ — ^§r- and the kinetic energy op- 
erator for the vacuum electron K = — ^- ^-if 2 (d~ + d~) . 

In the coordinates (x, y, 2) £ of @ becomes a plane 
with a constant x. We choose the separation point at 
the maximum of V(x), which is x p = . Thus it is 
easy to obtain using the JWKB approximation for 
the Green's functions. We assume the potential in the 
region beyond the anode to be a constant — V a . The 
emission wave arriving on the anode plane is the ana- 
logue of the two-slit Young's interference with an effec- 
tive slit constant z„ — < /— and an effective wave num- 



ber k a = j-y/2mV a . For a large focal length x/ and small 
emission angles, the wave function on the BFP is the 
Fourier transformation of the wave function on the an- 
ode plane that turns out to be 
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where 
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— and 6 Z = — are emission angles with 



respect to the y and z axes respetivelty, S p = (f^) 2 

and <P P {6 Z ) = ( 2 -^)hQ z e-^ e l with Mp = (1 + log(2 + 

V3))^k a Zp. The factors c ~ 0.890 and b = 1 + ^= for 
the present x p . Their values will change with x p but the 
effect will be compensated by the change of S p and the 
wave function is only slightly affected. 

With the emission wave function obtained, it is straight 
forward to find the emission current density of each 
individual graphene state. The total emission current 
is a non-coherent summation of the individual currents 
weighted by the Fermi-Dirac distribution(FDD). Because 
the emission depends strongly on the barrier height and 
is also controlled by the FDD, it is expected that only 
modes very close to the Fermi energy (Ep) can be emit- 
ted. This is the reason why CFE leads to observable 
interference even at high temperatures. Therefore the 
value of Ep is crucial for the emission current and should 
be determined taking account also of the induced ex- 
cess charge density. Using the density of states of the 
tight-binding model and an iteration method we find 

Ep = hvp^l 47V£ ° F "^ at the edge, consistent with a re- 



cent estimate. [13] The angular distribution of the emis- 
sion current on the BFP now reads 
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where E a = hvpk a . Because sin 7 = -j^0 y • - u y 
<t>p ~ Qzi it is clear that the current intensity is zero and 
the BFP image is dark in the center cross with y z = 
0. The pattern shown in FigfjTf) (and more in FigJ5]) 
is like a dragonfly. The dark slit that parallel to the 



graphene edge is caused by the destructive interference 
of two partial waves emitted from the p-orbital that has 
two lobes with the phase difference 7r on either side of 
the graphene plane (FigOTc)). Because (j> p only weakly 
depends on E(k) and 9 y , one may move it out off the 
integral ((3]) and see that the brightest point of a wing has 
an angle to the y-axis of fip 1 in which E p is substituted 
with Ep as the subscript F indicated. Clearly this angle 
is almost independent of the temperature. In Figj2l the 
image for T = 100 (a), 300K (b), and 1000K (c) are 
presented. In the low temperatures, the wing tips are 
sharply cut at 9 y ~ • In the high temperatures, 

the wing tips are smeared while the side-edges are still 
sharp. On the other hand, the dark "dragonfly" body 
has different originator. It is because that the angular 
distribution is proportional to sin 2 7 while 7, the mixing 
angle of the pseudo spin, is vanishing for the eigenstates 
with ki = 0. Before the lens where the waves of different 
&2 have not been dispersed, one will see two bright bands 
parallel to the y-axis(FigOJd)). 

The total emission current J may be given by the in- 
tegration of (U]). But it is more convenient to calculate 
it with the wave function on the anode plane. The J-Fq 
characteristic is presented in Fig |3] where the four curves 
from bottom to top are total currents for the temper- 
atures 100, 300, 600, and 900K, respectively. The cor- 
responding Fowler-Nordheim plots are also given as the 
inset of FigJ3J In the CFE regime one may expand the 
emission energy around Ep and obtain an explicit ex- 
pression for the total emission current, 
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0.092, a T = 5 ^w F and a F = 



-0 ~ - 2/3W F ailu "F ~ 8mv F W F with 

the subscript F indicating the substitution of E p by Ep 
in the corresponding quantities, and Io (Ii) is the modi- 
fied Bessel functions of the zero (first) order respectively. 
For the low temperatures and low fields such that both 
ax and apEp are small, the second line of (JSJ is not sen- 
sitive to either the temperatures or fields. The character- 
istic of the total emission current is mainly determined 
by the factor e~ b ° SF that is only slightly different from 
the conventional CFE theory. [23| Note that, as a charac- 
teristic feature of the two-dimensional CFE, J depends 
on h and Fq through Sf and Ep in a combination of 
F Vh. 

In summary, we have presented an analytical treatment 
of field emission from graphene at the microscopic level, 
including analytical expressions for the vacuum barrier, 
the Fermi energy, and the emission wave function in 
three-dimensional space. The Fermi energy is field de- 
pending due to the field penetration. The dragonfly-like 
emission pattern is predicted. The origin of the dragon- 
fly pattern is different from the interference of coherent 
beams from different atoms or atom clusters that give a 
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bright center, but is associated to the symmetry of the 
p-orbital of carbon and the bulk states of the graphene. 
The dark slit in the center of the pattern and parallel 
to the edge is a result of the parity of the p-orbital of 
the carbon atoms. We have proved that the width of 
the dark slit is almost independent of the temperatures. 
An less expected fact is the dark body of the "dragon- 
fly" that originates from the graphene structure. Intu- 
itively, the uniform edge of graphene plays the role of 
a one-dimensional electron grating. Along the edge of 
graphene, the transverse momentum of the emission is 
equal to the Bloch momentum of the graphene. The vac- 
uum barrier allows emission only from states within a 
narrow range of energy close to the Fermi energy and 
with small transverse kinetic energy. Because the mixing 
of electron states of the two sublattices, the emission of 
vanishing y-directional kinetic energy is absent, resulting 
the dark "dragonfly" body. Interestingly, the wing tips 
turn out to be brightest region. That implies the trans- 
verse kinetic energy is maximized for a given total kinetic 
energy controlled by Ep. The appearance paradox may 
be solved by the fact that the massless exciton of the 
graphene acquires a nonzero mass in vacuum. It is the 
massive real electron that tunnels through the vacuum 



barrier and its transverse kinetic energy is suppressed by 
the mass. Last but not least, the dark slit between a pair 
of wings on either side of the pattern may be smeared 
out by a magnetic flux parallel to the edge but the dark 
"dragonfly" body can not be changed. 

Our theory may get support from a recent experiment 
of Yamaguchi et al [25| who show two field emission pat- 
terns from a graphene edge in the configuration consid- 
ered by us. At a voltage of 5.1 V they see a regular 
interference pattern with the center the brightest that 
they attribute to interference between regular patches of 
oxygen atoms adsorbed on the edge. However, at a lower 
emission voltage of 2.1 V they see indeed a pattern that 
is dark in the forward direction flanked by two bright 
emission cones as predicted by our theory. The unusual 
up bending Fowler-Nordmeim plot (the inset of FigOJ) in 
the regime of weak fields and/or short graphene height 
has been confirmed by another experiment. [iH ] 
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FIG. 1: (a) A piece of graphene with the black balls repre- 
senting the carbons and the bond length a — 0.142 nm. The 
armchair edge is in the y-direction. The red rectangular pla- 
quette is a unit cell with the unit vectors ai = (V3a, 0) and 
a.2 = (0, 3a). (b) The energy bands of the Brillouin zone. The 
unit vectors of the reciprocal lattice are bi = (27r/(V3a), 0) 
and b2 = (0, 27r/(3a)). (c) The wave function profile of the 
HOMO along the armchair edge (DFT simulation), (d) The 
emission image near the anode plane in a field 7 ? o=20V//^m 
at the room temperature, (e) The electric objective lens, (f) 
The image on the back focal plane of the objective lens. 
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FIG. 2: (Color online) The image on the back focal plane for a 
graphene height h — 10/im and an anode distance x a = 1mm 
in the field Fo = 20 V / fim. The showing region has \9 y \ < 
0.005 and \9 Z \ < 0.00015. The panels (a), (b), and (c) are for 
temperatures 100K, 300K, and 1000K, respectively. 
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FIG. 3: (Color online) The total emission currents ver- 
sus the applied fields for h = 10/xm and an edge length 
3aN = 1mm. The curves from bottom to top are the to- 
tal currents for temperatures 100K, 300K, 600K, and 900K 
respectively. The inset shows the corresponding Fowler- 
Nordheim plots with ^/(^r) as the horizontal axis and 
log( JF^~ 2 )/(log(mA(nrn) 2 V~ 2 )) as the vertical axis. 



